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Abstract 

In this paper, we consider a unified approach to model wireless channels by a mixture of Gaussian 
(MoG) distributions. The proposed approach provides an accurate approximation for the envelope and 
the signal-to-noise ratio (SNR) distributions of wireless channels. Simulation results have shown that 
our model can accurately characterize multipath fading and composite fading channels. We utilize the 
well known expectation-maximization algorithm to estimate the parameters of the MoG distribution and 
further utilize the Bayesian information criterion to determine the number of mixture components auto¬ 
matically. We employ the Kullback-Leibler divergence and the mean square error criteria to demonstrate 
that our distribution provides both high accuracy and low computational complexity. Additionally, we 
provide closed-form expressions or approximations for several performance metrics used in wireless 
communication systems, including the moment generating function, the raw moments, the amount of 
fading, the outage probability, the average channel capacity, and the probability of energy detection for 
cognitive radio. Numerical Analysis and Monte-Carlo simulation results are presented to corroborate 
the analytical results. 
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I. Introduction 

Modeling the terrestrial wireless propagation is of importance for the design and performance 
analysis of wireless systems. In a typical mobile radio propagation scenario, the received signal 
presents small scale power fluctuations, due to multipath propagation, superimposed on large 
scale signal power fluctuations, also known as shadowing, which is due to the presence of 
large obstacles between the transmitter and receiver. The small scale fading results in very 
rapid fluctuations around the mean signal level, while shadowing gives rise to relatively slow 
variations of the mean signal level [1]. A common example of composite fading channels is the 
Nakagami-Lognormal (NL) channel. In this case, the density function is obtained by averaging 
the instantaneous Nakagami-m fading average power over the conditional probability density 
function (pdf) of the log-normal shadowing, resulting in a complicated pdf that has no closed 
form expression [2]. 

The K [3] and generalized-/! (Kg) distributions [4], [5], have been introduced as relatively 
simpler models to characterize composite fading channels, in which the Lognormal distribution 
is replaced by the Gamma distribution in the Rayleigh-Lognormal (RL) and NL distributions, 
respectively. The k — /i and the // — // distributions [6] are general fading distributions for line- 
of-sight (LOS) and non-line-of-sight applications, respectively. These distributions can represent 
the Rice (Nakagami-n), the Nakagami-m, the Rayleigh, the One-Sided Gaussian, and the Hoyt 
(Nakagami-g) distributions as special cases. Quite recently, a generalization of the n-fi fading 
channel, where the LOS component is shadowed, has been proposed in [7]. All these models 
contain the modified Bessel function of the first or second kind, which complicates further 
analytical performance measures. In [8], the Lognormal distribution was replaced by the Inverse- 
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Gaussian distribution, resulting in the Rayleigh/Inverse Gaussian (RIGD) distribution, followed 
by its generalized versions, i.e., the (^-distribution [9], the k — ^/Inverse Gaussian distribution 
[10] and the r] — ^/Inverse Gaussian distribution [11]. The drawback of these distributions is 
their increased complexity due to the presence of the modified Bessel function of the second 
kind. Recently, an interesting work has been proposed by Atapattu et al. [12], where several 
channel models were expressed as a mixture Gamma {MQ) distribution via Gauss-Quadrature 
approximation. The MQ model is more accurate than the aforementioned alternatives, and it has 
the advantage of simplicity as well. 

Finite mixtures of distributions provide a mathematical-based approach to statistical modeling 
of a wide variety of random phenomena [13]. In this paper, an alternative model, that represents 
both composite and non-composite fading channels by Mixture of Gaussian (MoG) distributions 
is presented. The approximation method is based on the expectation-maximization (EM) algo¬ 
rithm, which was coined by Dempster et al. in their seminal paper [14]. The EM algorithm is 
essentially a set of algorithms exceptionally useful for finding the maximum likelihood estimator 
(MLE) of any distribution in the exponential family [15], and widely used for the missing data 
problem (i.e., modeling a mixture distribution). The main contributions of this paper can be 
summarized as: 

• We propose MoG distributions to model both the envelope and the signal-to-noise ratio 
(SNR) distributions of wireless channels. The proposed approach is proven to accurately 
model both composite and non-composite channels in a very simple expression. 

• We determine the number of components using the Bayesian information criterion (BIC), 
while the corresponding parameters for the mixture are evaluated using the EM algorithm. 

• We demonstrate the importance and tractability of our model by deriving several tools for 
the performance analysis of single-user communications such as the outage probability and 
raw moments. Moreover, we derive the moment generating function (MGF), of which the 
symbol error rate (SER) of L-branch maximal ratio combining (MRC) diversity system is 
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presented for various signaling schemes. Finally, we derive an approximation for the average 
detection probability in cognitive radio networks and provide an upper bound to the error. 

• Numerical analysis and Monte Carlo simulation results are presented to corroborate the 
derived analytical results. 

The rest of this paper is organized as follows. Section II gives a brief description of several 
wireless channel models of interest. In Section III, the MoG distribution is introduced together 
with a brief description of the EM algorithm. Section IV presents a detailed comparison of 
the MoG distribution to the channel models it can approximate. In Section V, performance 
metrics, such as the MGF, the raw moments, the amount of fading, the outage probability, and 
the average channel capacity are derived using the MoG distribution. Simulation results and 
numerical analysis are presented in Section VI. Finally Section VII concludes this work. 


II. Fading Channels 


Radio-wave propagation through wireless channels undergoes detrimental effects characterized 
by multipath fading and shadowing. Modeling of such fading channels is typically a complex 
process and often leads to intractable solutions. Considerable efforts have focused on the statis¬ 
tical modeling which resulted in a wide range of statistical models for fading channels [2]. 

In this section, we give a brief description of some well known channels that often lead to 
intractable performance analysis of wireless communication systems. 

A. The Nakagami-Lognormal Channel 

The NF fading model is a mixture of Nakagami-m distribution and Fognormal distribution 
obtained by averaging the instantaneous Nakagami-m fading average power over the pdf of the 
log-normal shadowing as follows 



( 1 ) 
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where / a (a|cr) is the Nakagami-m distribution given by 


f Z L 2m— 1 —m— 

Ta(ct\cr) = ——rOt e 

J y 1 ' cr m T (m) 


( 2 ) 


where T(.) is the gamma function [16] and m is the fading parameter, which is inversely 
proportional to multipath fading severity i.e., as m —» oo, multipath severity diminishes. Note 
here that m = ^',7, . The average power a follows a Lognormal distribution, contributing to 
shadowing at longer routes, expressed as 


fo- (o’) = 


-(101og(cr )-M) z 
l 2C 2 


( 3 ) 


\f2Hcr\C, 

where A = M and ( 2 , measured in dB, are the mean and variance of the Gaussian RV 
V — 10 log 10 (cr), respectively. In order to compare (2) with that of the Gaussian RV X = ln(cr), 
the following relations apply [17] 

X = XV, (4) 

M x = AM, 

Cx = AC- 

An important remark regarding the Lognormal distribution is that, while C essentially defines 

M 

different Lognormal distributions, M is effectively a scaling factor [17]. Denoting M n = 10 io 
and x — jj-, then it is straightforward to show that 


fa ( a M n ) — —/„ (a\M — 0). 


( 5 ) 


Therefore, it is only sufficient to perform an approximation for M = 0 dB, and generalize the 
results to other scaling factors. 

Let E s denote the energy per symbol, N 0 be the single sided power spectral density of the 
complex additive white Gaussian noise (AWGN) and assuming E [|cr:] = 1, where E[.] denotes 
the expectation operator. By applying the following transformation to (1) 


7 


= « 2 7> 


( 6 ) 
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where 7 = £[ 7 ] = ^ is the average SNR, we obtain the NL distribution as 


fi 0 ) = 


m no o m—1 

/ X 


(8.686) m' 

T (m) \/2vrC Jo 7 




e ^ e 2 c 


da. 


(7) 


The SNR density function is not expressed in a closed form, making the performance analysis 
of wireless communications under this particular channel very complicated or intractable. Note 
that the RL distribution is a special case of NL distribution with m = 1. 


B. The k — /i and 77 — /i Fading Models 

With the emergence of relatively new wireless mediums, such as which occurs in underwater 
acoustic [18], [19] and body communication [20]-[23] fading channels, the wireless research 
community have had a reincarnated interest in finding more accurate and generalized fading 
models that provide a better fit to new and realistic measurements. Consequently, new fading 
models, such as the n — /i and 7 — 77 distributions, were proposed [ 6 ]. The n — /i fading model 
is mostly used to represent the multipath fading with LOS condition and includes the following 
fading models as special cases: the Nakagami-n (Rice), the Nakagami-m, the Rayleigh, and the 
One-sided Gaussian. The instantaneous k — /i SNR distribution is expressed as [ 6 ]. 


/ \ M+1 

Cl+rA 2 


n -7 


fi (a) = -g=r- 7—7 

n 2 exp {/in) 


m-i f /i (1 T k) 


x 2 exp 


—- X ) 4_i I 2 fly 


' K (1 + K) X 

7 


( 8 ) 


where At > 0 is the ratio between the total power of the dominant components and the total 
power of the scattered waves, /1 > 0 is given by /a = V | ', 7 ^ | ^ 2 . and I fl (.) is the modified 
Bessel function of the first kind and order /i [16, eq. 8.445]. It is worth mentioning that as k 
tends to zero, the k — /i distribution degenerates to the exact Nakagami-m distribution, with 
/i = m = 77777 - Additionally, by setting fi = 1, the k — fi distribution degenerates to the exact 
Nakagami-n distribution, with n = n. 

Complementing the k — // model, the // — // model was proposed to represent NLOS multipath 
environments, where it includes the Nakagami-g, the Nakagami-m, and the One-Sided Gaussian 
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distributions as special cases. The instantaneous 77 — n SNR distribution is expressed as [6]. 


2- v /7T/i At+ 2 

4 \ X ) = -7-— X' 




exp 


2 fihx 


4-i 


/ 2 jiH 


-x 


( 9 ) 


r (/i) ‘ V 7 / ^ 2 V 7 

where /1 > 0 is given by n = |+ 7; 77 } (-f-) 2 , and parameters h and H can have two 
different formats corresponding to two different physical phenomena as follows: In Format 1 , 


h = 


2+r ? ~ 1 +7 ? 


and H = 


v -v 


, where 0 < 77 < 00 is interpreted as the power ratio between 


the independent in-phase and quadrature components. In Format 2 , the in-phase and quadrature 
components are correlated and have a power ratio of unity. The two corresponding parameters 
are defined by h = j ^4 and II = where — 1 < 77 < 1 represents the correlation between 
the in-phase and quadrature components. The two formats can be obtained from each other using 
the relation T/pormati = ] • 0' s worth mentioning that the 77-/7 distribution degenerates 

to the Nakagami-g distribution by setting /i — 0 , with q = in Format 1 and q = !—- 

Format 2 . 


l+r) 


in 


C. The k — ji Shadowed Fading Models 


The k — n Shadowed fading model, was firstly proposed as a LOS Shadow fading model 
[ 7 ], where unlike the NL formulation above, it is assumed that only the dominant components 
of the multipath clusters are subject to random fluctuations. The unconditional instantaneous 
SNR distribution of the ft — /1 Shadowed model is obtained by averaging the conditional n — // 


distribution over the Nakagami-m distribution as follows 

/OO 

41« (7; 0 4 (0 4 




/i (1 + k) 2 

’ Z M-i 

7 k 2 


M-l 
\ 2 

7 


exp 


/i( 1 + k) , j\ m 


7 


( 10 ) 


0 ( 7 ), 


where 0 (7) = f 0 °° 2 exp (-£ 2 (/xk + m)) £ 2m 1 (2^ 

following closed-form expression 


k( 1 +b )7 
7 


r (m) 

d£, which results in the 


4 (7) = 


(1 + /7 

F {fJt) 7 + m) m 77 


M-l 


exp 


/4 1 + 47 

7 


1-4 m, fi\ 


/7 2 7C(1 + K) 7 
fin + m 7 


, ( 11 ) 
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where the function ] F ] (.,.) is the confluent hypergeometric function [16, eq. 9.210.1], defined 

nc F (n h- ~) — 1 4- a ( a +l) z 2 , a(a+l)(a+2) z 3 , 

V u > U ’ *9 1 6 1! ” r 6(6+1) 2! ^ 6(6+1)(6+2) 3! ' "" 

Interestingly, in a very recent work [24], it has been shown that under a new formulation the 


K — fjL Shadowed fading model can represent the // — // distribution as a special case with /i = 2/z, 
k = and m = m, where the underlined symbols belong to the k — // Shadowed model for 
the sake of clarity. 


III. The MoG Distribution 

We consider the problem of estimating the wireless channels’ density functions. Gaussian 
mixtures [25]-[29] are often used due to the fact that their individual densities are efficiently 
characterized by the first two moments [30], [31]. The MoG distribution is attributed to have 
the Universal-approximation property, as it has been proven by Weiners approximation theorem 
[25], which states that the MoG distribution can approximate any arbitrarily shaped non-Gaussian 
density. The objective of this section is to provide a unified MoG distribution that can accurately 
represent different fading channels. 

A. Parameter Estimation of the MoG Distribution 

Let the i th entry of a random data vector Y = (yi,..,y n ), which represents the channel fading 
amplitude of the composite models, be regarded as incomplete data and modeled as a finite 
mixture of Gaussians as follows 

c 

p(Vi\9) = ( 12 ) 

3= 1 

where i = 1,..., n and C represents the number of components. Each j th component is expressed 
as 

<13) 
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where the weight of the j th component is ujj > 0, with Uj = 1. The parameter 6j = ( fi 3 . rjj j 
correspond to the mean and variance of the j th component, respectively. 

Let the complete data X be the joint probability between Y and Z, where Z e {1, ..,C} is 
a hidden (latent) discrete RV that defines which Gaussian component the data vector Y comes 
from, namely, 

P(Z = j) = j = (14) 

Ideally, one would like to maximize the log-likelihood function as follows 


Omle 


argmaxL (6) 
ee© 


(15) 


= argmaxlog p (y\0). 


However, maximizing L (0) is not tractable and difficult to optimize [32]. Instead, the EM 


algorithm solves the MLE problem by maximizing the so-called Q-function as follows [33] 


9 (m+i) = argmaxQ (0|0M) 

= arg max ’E x \ v fl^) [ lo S Vx {X\9 )], 


(16) 


where m is the iteration index. The EM algorithm is performed by two iterative steps, namely 
the expectation step (E-step), and the maximization step (M-step). We set initial guesses of the 
MoG coefficients, i.e. a/°\ /i®, r/°\ whereby in the E-step, we compute the posterior probability 
(membership probability) 


(m) 

Pa = 




(17) 


In the M-step, the coefficients are updated by differentiating the Q-function with respect to ui, 
/i, and r], resulting in the following analytical (m + 1 ) lh estimates 


to 


,, (m+l) 


(m+1) _ 1 (m) • _ 


n 




1=1 

n 


j = 


n 


(in) 


(18) 

(19) 


i=l 
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( 20 ) 


This iterative procedure is terminated upon convergence, that is when [J rn + V > — £( m )| < S, 
where 



( 21 ) 


is the log-likelihood, and <5 is a preset threshold. 

The EM algorithm is guaranteed not to get worse as it iterates, i.e. < L ( m ) [14]. 

Hence, the lower 5 is set, the more accurate the approximation would be. In addition, one can 
always increase the accuracy by increasing the number of components. Though, this technique 
might be stuck in a local maxima, since the likelihood is a marginal distribution. However, one 
could mitigate this problem by heuristics and multiple initial guesses. In this regard, Do et al. 
[32] suggest to initialize parameters in a way that breaks symmetry in mixture models. Finally, 
it is noteworthy to point out that the EM algorithm has an advantage of being a completely 
unsupervised learning algorithm, which makes it very convenient for our density estimation 
application. For more details, one can refer to [15], [33], and references therein. 

B. The pdf of the Instantaneous SNR of the MoG Model 

With the aid of the EM algorithm, all fading channels’ amplitudes can be represented as 



( 22 ) 


By taking the change of variables 7 = 7X 2 , the pdf of the instantaneous SNR of the MoG 
distribution can be written as 
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C. Determining the optimal number of mixture components 

Generally, when fitting a finite mixture distribution, the determination of an appropriate number 
of mixture components is inevitably a necessity. Choosing a small number of components would 
yield an inaccurate representation, while a very large number of components would unnecessarily 
increase the complexity of the distribution and may cause over-fitting. In addition, Chen [34] has 
shown that knowledge of the the number of components yields a faster optimal convergence rate 
for the estimates of a finite mixture than it would when the number of components is unknown. 

In this subsection, in order to derive an appropriate number of mixture components, we adopt 
a simple yet effective unsupervised information theoretic criterion, called the BIC, which was 
introduced by Gideon Schwarz in [35]. 

Let x = {xi, ...,x z , ..xm}, correspond to M independent and identically distributed ( i.i.d.) 
samples, drawn from any of the envelope distributions of the actual aforementioned fading 
models, then the log-likelihood function of the MoG distribution can be expressed as 

M ( C 


L c (fl)=lnPr(*|e,C) = ^;in]E v/ 2i j)? 


Uli 


exp 


2 = 1 


[Xz ~ hi) 
2f] 2 


(24) 


where 0 = [(£>1,0)2, ...,(£><7, A ■■■,frc,V 1 , 7 ) 2 , ■■■,Vc\ are the estimated parameters and C is the 
corresponding number of components. The corresponding BIC score can be computed as 


BIC c = -2L r (e 


c 


C ln(Af). 


(25) 


It can be seen that the BIC penalizes the model complexity by adding the regularization 
coefficient, C In (At). It is worth noting that although the EM algorithm maximizes the log- 
likelihood distribution, the BIC is an asymptotic approximation to the transformation of the 
Bayesian a posteriori probability, Yr(6\x 1 C). As such, in a large-sample setting, the number 
of components determined by the BIC is asymptotically optimal from the perspective of the 
Bayesian posterior probability. Here, we select the candidate model satisfying the minimum BIC 
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score, satisfying asymptotically the maximum Bayesian posterior probability as 

C opt = arg min BIC C - (26) 

Ce N 

Fig. 1 depicts the normalized BIC versus the number of components for some fading scenarios 
selected from Section II. The corresponding optimal number of components, C op t, indicated in 
the legend, will be adopted in the simulations and numerical results hereafter and will be denoted 
by C. 



Figure 1. Normalized BIC versus the number of components 


Fig. 2 shows the optimal number of components as a function of the amount of fading. It 


is observed that as the fading becomes more severe, the mixture requires more components to 


accurately characterize the channel. 
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Figure 2. Optimal number of components versus the amount of fading 


IV. MoG Model Analysis and Comparisons 

In this section, several scenarios of the aforementioned fading channels are approximated using 
the MoG distribution, as in (23). The number of components was selected automatically using 
the BIC method explained in Section III-C. We point out that higher accuracy can be achieved by 
increasing the number of components. In order to validate the accuracy of the approximations, 
we use two criteria of error, namely the MSE and the KL, defined as 


MSE 


E 


AO) 



(27) 


and 

KL ff 7 (x), f 7 (x )) = j / 7 (x) log dx, (28) 

7 o J7 

respectively. Here, / 7 (x) is the exact pdf, and / 7 (x) is the approximated pdf (MoG). The KL 
divergence, also known as relative entropy, is an information theoretic measure that quantifies 
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the information lost when / 7 (x) is used to approximate / 7 (x) [36]. Note that the MSE and KL 
measures are used in several related works, see e.g., [ 8 ], [9], [12]. 

Fig. 3 provides the approximation results for several scenarios of the NL, k — /i, 77 — /i, and 
k — fi Shadowed fading models. The corresponding number of components are indicated in the 
legend. As shown from the MSE and KL measures, the approximation is very accurate when 
both increasing the shadowing and the multipath fading severity, whereas for large amount of 
fading, the number of components increases. For instance, for the 77-/1 distribution, when 77 = 5, 
and /i = 10, the amount of fading and number of components were 0.0721 and 3, respectively, 
whereas when 77 = 0.7 and // = 0.4, the amount of fading and number of components increased 
to 1.2999 and 8 , respectively. The parameters of the approximations are tabulated in Appendix 
A. 



O m=1, C=3dB, C=7, KL=2.2e-4, MSE=8.0e-6 

* m=4, C=1dB, C=4, KL=1.1e-4, MSE=8.3e-6 
O k=1 , p=1, C=6, KL=1,6e-4, MSE=1,3e-5 

* k=3, p=1, C=5, KL=1,2e-4, MSE=1,0e-5 

O ri=0.7, p=0.4, C=8, KL=6.1 e-4, MSE=6.9e-5 

* il=5, p=10, C=3, KL=1.1e-4, MSE=2.3e-5 
k=1, n=3, m=3, C=4, KL=1.2e-4, MSE=9.8e-6 

-Exact Distribution 


Figure 3. MoG approximation for different channel models. 


Further verification of the accuracy via numerical means is addressed in section VI. The 
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purpose of this approach is not to increase the accuracy of the approximation, as it has already 
been achieved in all aforementioned fading alternatives, but rather to provide another unifying 
and simplifying distribution that has the potential of approximating all contemporary fading, 
composite and non-composite, models via the EM algorithm. 

V. Performance Analysis of Wireless Channels 

The MoG distribution provides a simplifying and unifying analysis for wireless communi¬ 
cation systems over various composite and multipath fading channel models. In this section, 
we first derive several performance metrics, which can be used for the evaluation of wireless 
communication systems in a generalized manner. In particular, we derive expressions for the 
raw moments of the MoG model, the amount of fading (AF), the outage probability, the channel 
capacity, and the MGF. We further derive expressions for the SER performance of L-branch 
MRC diversity system and the average probability of detection for cognitive radio networks. 


A. Moment Generating Function 

By definition, the MGF M 7 (s) = E [e~ sl ] is given by 


c 


i=l 


r ° 1 

( ( 


/ exp 

Jo Vl 


2 V- 

/ 


e 1S d'y. 


(29) 


Applying the change of variables x — ,/=, and after expanding the exponentials and considerable 


mathematical simplifications, we get 

c 


M M = Ey=^/ o ex P | — — ]dXt 


2 Vt 


where = 1 + ‘hj-xs. Then, after some mathematical manipulations, we obtain 


(30) 


M, (s) = Y, 




1=1 


VAvr 


exp (— z 2 ) dz, 


(31) 
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which leads to the following expression 

c 


M y (s) = V —4= exp 

7 trVA 




Q 


Hi 


Pi ) V ViVPi) ’ 

where Q(.) is the Gaussian Q-function defined as Q (x) = -P= J x °° exp 


-f)du. 


(32) 


B. Raw Moments 


The n m raw moment of the MoG distribution, by definition, is 

( ' ^ ' 2 


c rc 

e [ 7 n ] = y Jg_ / 
V&Ynm J 0 


r-oo 

VP 


exp 




V 




d'y. 


(33) 




By taking the change of variables x — ./=, and after some mathematical simplifications, we get 


c 

em = I>r / 

i=i 7o 


°o ^, 2 n 


Alternatively, we can write (33) as 


\/27T rji 


c 


(x - Hi ) 2 

6XP 1 -V— I dx ’ 


(34) 


E [7 n ] = E [Af*] , 


(35) 




where Xi ~ N {mXi) is the Gaussian component. Using the MGF approach, (35) can be 
expressed as 

JL ^( 2n ) Mv r.«t) 

(36) 


EMI = E^M-7T75§—1.=0, 


Z =1 


ds( 2 ”) 


where Mx t (s) = E{e is the MGF of Xi, which is given by 


M x . (s) = exp /i;s + 


v y 


(37) 


Equation (36) is mathematically convenient for solving the first few moments. 

An alternative approach that yields a closed form expression can be attained by following the 
same method in [37], where the Ah raw moments of X, are derived as 

v 1 Vi 


r (- + -) 

E[X?] = rj?2* [2 V 2 iFi 


7r 


2 ’ 2 ’ 277 ?J ’ 


(38) 
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where v is an even integer (note that there is no loss in generality). 

By substituting (38) into (35), the n th raw moment of the MoG distribution is derived as 


c 


E [ 7 "] = 1F1 


i =1 


7r 


pi 


2 ’ H\ ’ 


(39) 


C. Amount of Fading 

The AF measure was firstly introduced by Charash [38], as a measure of the severity of 
the fading channel. The AF requires the knowledge of only the first two moments in the 
corresponding fading channel, where it is defined by 

E[7 2 ]—E[ 7 ] 2 


AF = , 

EM 

By solving (39) for the first two moments, we obtain 


(40) 


AF = Elii Ui (pf + 6/ifr/j 2 + 3 rjf) _ 1 




( 41 ) 


D. Outage Probability 


The outage probability is a standard performance criterion used over fading channels. It is 
defined as F(y th ) = fj"' f, (x) dx. By performing the following change of variables applied to 
(23) 


y = 


= - Pi 


Vi 


(42) 


and after some mathematical manipulations, the CDF of (23) can be written as 

C „_v £_ , O' 

r 

_ 2 

*— 1 v m ' 

Further simplifications yield 


Vi 


(43) 


c 

F ( 7 th) = Ui 
1=1 




1 1th _ ,, 

Vi 


(44) 
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E. Average Ergodic Channel Capacity 


When only the receiver has knowledge about the channel state information (CSI), the ergodic 
capacity C erg is expressed as 


jd roo 

Cerg = / In (! + 7) /t (7) d'y. (45) 

where B is the channel bandwidth measured in Hertz. Unfortunately, the exact solution of (45) 
is intractable. Instead, a computationally simple and very accurate form can be obtained by 
following [39], where ln(l + 7 ) is expanded about the mean value of the instantaneous SNR, 
£[7], using Taylor’s series, yielding 


-(1 +T ) = ln(l + K[ T ]) + E ( T fcgS 

Ul = l v L 


In (1 + E [7]) + E j 7 lf + o [(x — E [ 7 ]f] 


(46) 


1 + E [ 7 ] ' 2 (1 + E [ 7 ]) 2 
where o (.) is one of the Landau symbols defined as / = o (</>) means that ^ —> 0. 

Substituting the logarithm function approximation from (46) into (45), a second order approx¬ 
imation for the ergodic capacity can be evaluated. It is noted that the solution to (45) is obtained 
by taking the expectation of In (1 + 7 ). Hence, taking the expectation of (46), we get 


67 

er9 ~ In 2 


In (1 + E [ 7 ]) 


where E [ 7 "] is evaluated using (39). 


E[7 2 ]-E 2 [ 7 r 
2 (1 + E [ 7 ]) 2 J 


(47) 


F. Symbol Error Analysis 

In order to further demonstrate the significance of the MoG distribution, we study the perfor¬ 
mance of independent but not identically distributed ( i.n.i.d .) L-branch MRC diversity receiver 
over various composite and non-composite fading scenarios. The MRC scheme is the optimal 
combining scheme at the expense of increased complexity, where the receiver requires knowledge 
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of all channel fading parameters [2]. Here, the receiver sums up all received instantaneous SNR 
replicas as follows 

L 


2 MRC 



(48) 


k =1 


The corresponding MGF is thus 

L 

M W M = El- ^ } = If (s), (49) 

k =1 


where M 7fc (s) is derived in (32). The SER, P s (E), for coherent binary signals, can be computed 
as follows [2] 


P.(E) 


E 


Tmhc 


Q 



(50) 


where g is some constant resembling several coherent binary signals, such as coherent binary 
phase shift keying (BPSK) and coherent orthogonal binary frequency shift keying (BFSK) 
corresponding to g = 1 and g — \ , respectively. By substituting the Q-function by its definition 
in [2, eq. 4.2], the SER is written as 


P. ( E ) = - 

7r 


exp 


'o jo 


9 7 a 


sin 2 ( 9 ) 


f'YMRC (7 MRC )d'f A 


dO. 


The inner infinite integral in (51) is the equivalent MGF derived in (49), yielding 


P, (E) = - 

7T 


, ZL L 
2 


n w 


7fc 


k =1 


9 


sin 2 (6) 


dd , 


(51) 


(52) 


where M 7fe (.) was derived in (32). Following a similar approach, and by utilizing [2, eq. 8.23] 
and [2, eq. 8.12], the SER expressions for M-PSK and square M -QAM signaling schemes are 
given, respectively, by 


Ps(E) 



( sin2 (§) 
^ sin 2 ( 9) 


d6, 


(53) 


Ps(E) 


4 ( yfM- 1 
vr \ VM 





fc=i 


/ 9qam \ 
Vsin 2 (0)J 




>/M- l\ 




n v ' 




fe=l 


/ 9qam \ 
Vsin 2 (0)J 




(54) 


? 


where 9qam = | ( M - !)• 
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G. Probability of detection 


Cognitive radio (CR) is a promising technology that can enhance the performance of wireless 
communications [40]. The basic concept behind opportunistic CR, is that a secondary user 
is allowed to use the spectrum, which is assigned to a licensed primary user (PU), when 
the channel is idle [41]. The CR users perform spectrum sensing in order to identify idle 
spectrum. Energy detection is the most common sensing technique in CR networks, due to its 
low implementation complexity and no requirements for knowledge of the signal [42]. Several 
studies have been devoted to the analysis of the performance of energy detection-based spectrum 
sensing for different communication and fading scenarios [43]. The probability of detection, P d 
is an important performance evaluation metric representing the probability that an active PU is 
detected by the CR node. Over fading channels, the average detection probability is evaluated 
by averaging the AWGN channel P d over the SNR distribution as follows: 


Pd = j Qu (v^, V\) fry (7) d'y, (55) 

0 

where Q u (\f2r<, a/A^ is the probability of detection over AWGN channel with Q u (.,.) being 
the generalized Marcum-Q function [44]. Moreover, A is a predefined energy detection threshold, 
u is the time bandwidth product which corresponds to the number of samples of either the in- 
phase (/) or the quadrature (0) component, and 7 = AA- is the received SNR, where E s is the 
signal energy, N 0 the one-sided noise power spectral density, and the channel gain E [|a 2 |] = 1. 
The Marcum Q-function can be expressed as [45] 


Q u (^,V\)=e-’J 2 


T "r(« + r»,f) 


—' n\ T (u + n) 

—n \ / 


(56) 


n =0 


By substituting (23) and (56) in (55), we obtain 

w- e ‘-it 0 4 tA+ 2 ) T (u + n, |) T (2 n + 1) (777?) 


2—1 


\/2vr7 r\i 


E 


2n+l 

2 


2n+l 


-D 


=0 n!T (u + n) + 1) 2 


-(2n+l) 




jii\/ 2 itfi +1 

(57) 
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where D n (i) is the parabolic cylinder function [16]. 

The truncation of the infinite series in (57) generates the following error 

r(M+n,|)r(2n+l)(7r ?I 2 )~2— 


, _ Wi 

e t - Z-/»=i — 


Vi 

2 Vi - 


4~tvf +2 




Eoo 

n— 


oo 
n=p+l 


n!r(ii+n)( 27 ? 7 ?+l) 


Z-^i=l 


-JH 

2r h‘ 2 


4-17 ?+2 


V 2n 'IVi 


271 + 1 


2n+l 


-(2n+l) 


v—ioo r(u+n, 2 )r( 2 n+l)( 7 ry 2 ) 

2_/r*=0 77 i ■ 2 Ti+i r^-(2n+l) 


n!r(«+n)( 27 r;?+l)' 


2n + l 


r(**+n,f)r( 2 n+l)( 7 *j?)' - 

*£**=0 77 i n . 2 ra+i (2n+l) 


Pi 


Pi 


Pi\/27?7?+T 


(58) 


Pi 


riiy/2^+1 


Vi\/2 IVi+l 


(59) 


n!r(«+n)( 27 T; 2 +l) 

Since the T (a, x) function is monotonically decreasing with respect to x, thus, T (u + n, < 
T (u + n). Hence, the infinite sum in (59) can be upper bounded by: 


r = 


£ 

71=0 


T(2ra + 1) 


n\ 


W 


27 Vt + 1 


2n+l 

2 


L>_ 


(2n+l) 




ViV 2 W + 1 , 


(60) 


By expressing the parabolic cylinder function in terms of the confluent hypergeometric function 
according to [16, eq. (9.240)], one obtains 


4t7(2-t7 i 2 + 1 ) 


t = yfire '■*' '■'* 


T(2n + 1 ) ( yrj, 

/ v lo 2n + 1 

n=0 n!2 2 


9 \ 2 n±i 

2 \ 2 


2 77* 2 + 1 
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n\ 


11 

i Fi [n + —, —, 






i F\ ( n + 1 , —, 


7+ 


2’ 2’ 2^(2^ + 1)7 T (n + 1) vW 2 W + 1 V ’ 2 ’ W( 2 W + 1) 

( 61 ) 

By expanding the involved hypergeometric functions and after some algebraic manipulations, 
the following equality is valid 


T =y /^ e ~4^A+D ^—p (n + 2)T(n + 2) 


n =0 2=0 
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477* + 2 
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(62) 
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By recalling that x = x\ j(x — 1)! and with the aid of the Pochhammer symbol identities, it 
follows that 


T =1 


'wm 


47 Vi + 2 


477 ^( 2 — 77 ^+ 1 ) 


iv 2 


(Dn+i fi)n ( 4yn'i+ 2 ) \ 2 Vi( 2 r r ii+ 1 ) 


SS GJiWn n! 
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t 1 


V2 


iiy/ 2r rVi + 1 l^o V 2 


EEtS 


7+ 2 


(1 )n+i (^)n 1 477 ^+ 2 ) l 277 ^( 27 ^ + !) 


)i(^ 


n! 


11 


(63) 


Importantly, the above infinite series can be expressed in closed-form in terms of the Humbert 
'Pi function [46] yielding 
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e* < E 


Wie 2r n 




K 


^ yJl 3 - 1 -jibL _ 

+ 2 1 \2’ ’2’ ’477/1+ 2’2^(27^ + !) 
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-£> 


=0 7 i!r (w + n) (27772 + 1 ) 2 


-(2n+l) 




7/^27772 + 1 


(64) 


VI. Simulation Results 

In this section, we present some analytical and simulation results for the outage probability, 
the average ergodic capacity, the SER of MRC scheme and the average detection probability for 
cognitive radio. 

Fig. 4 and 5 depict the outage probability, as in (44), versus the threshold SNR ^ th and the 
average SNR 7, respectively. Two NL scenarios are considered in Fig. 4, where the multipath 
severity is reduced from 777 = 1 to m = 3 , and two k — n scenarios are considered in Fig. 
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5, where the control parameter k is increased from k — 1 to k — 3. Here, one can notice 
how accurate the approximation is. Also, it is very noticeable how the multipath fading severity 
affects the outage probability performance. 



8 10 12 
Threshold SNR y (h (dB) 


(Analytical) m=1,£=1 dB 
+ (Analytical) m=3,£=1 dB 
Simulation 


7 


X 

= 20dB 


Figure 4 . Analytical and simulated outage probability versus 7 th for two scenarios. 


Fig. 6 depicts the capacity, as in (47), for some selected scenarios from the 77 — /x, k — // 
and k — /i Shadowed fading models. As shown, the severe NLOS configuration of the r/ — // 
scenario exhibits the worse capacity. In addition, one can see that introducing the Nakagami-m 
shadowing to the k — /1 distribution has worsened the capacity. The term simulation refers to 
cross-validating the results via the recursive adaptive simpson quadrature method performed by 
the aid of a mathematical package. 
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Figure 5. Analytical and simulated outage probability versus 7 for two scenarios. 



Figure 6. Analytical and simulated ergodic capacity, B=^. 
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Fig. 7 illustrates the analytical SER of the BPSK signaling scheme for various NL scenarios, 
including mild and severe fading cases in both shadowing and multipath. The solid squared line 
represents the corresponding Monte Carlo simulation. It is quite noticeable how the multipath 
severity plays greater role in determining the SER, where as observed, incrementing m by only 
1 yields a SER performance improvement of about an order of magnitude at observed mid-range 
average SNR values. On the other hand, increasing ( from 1 to 3 dB, while fixing m, yields a 
very similar SER performance. 



Average SNR 7 (dB) 


(Analytical) m=1,£=1 dB 
(Analytical) m=1 .£=2 dB 1 • 

(Analytical) m=1 .^=3 dB 
(Analytical) m=2.i)=1 dB 
(Analytical) m=3.i(=1 dB 
Monte Carlo Simulation 


Figure 7. Analytical and simulation SER of 2-branch MRC diversity receiver for BPSK signaling scheme for RL and NL 
fading channels. 


Fig. 8 features the analytical SER of the 16-QAM signaling scheme for various k—/i Shadowed 
fading scenarios, where the control parameter /i and the shadowing severity parameter m are 
varied. Here, L corresponds to the number of antennas, and it is noticeable that our model is 
still very accurate for high antenna diversity order. 
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(Analytical) K=1,n=3,m=1 (L=1) 
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0 

(Analytical) k=1 ,h= 1 ,m=3 (L=2) 

-B— 
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15 


20 


Average SNR y Q (dB) 


Figure 8 . Analytical and simulation SER of L-branch MRC diversity receiver for 16-QAM signaling scheme for for various 
k — fi Shadowed fading scenarios. 



Figure 9. ROC curves for Nakagami-Lognormal channel with u = 3 and 7 = 5 dB 
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Fig. 9 depicts the receiver operating characteristic (ROC) curve, where P d evaluated using 
(57) is compared to the theoretical expression from (55) substituting / 7 (7) by the NL fading 
channel pdf. These curves, depicting P d versus the probability of false alarm (Pj), are essential 
for the performance evaluation of the CR. We consider a single user performing spectrum sensing 
over NL channel where m is the Nakagami-m fading parameter and ( 2 , measured in dB, is the 
variance of the Gaussian random variable defined by V = lOlogio(a), where a corresponds 
to the Lognormal shadowing. The parameters for the MoG distribution are calculated via the 
EM algorithm and are provided in Appendix A. It can be observed that, for both cases, the 
approximation is very accurate over all the values of {Pj). Moreover, by evaluating the truncation 
error upper bound using (64), we conclude that although the truncation error depends on many 
parameters, considering the practical case where Pf = 0.1 and truncating the series at n = 12 
ensures that the error does not exceed 0.08 for the considered channels. 

VII. Conclusion 

In this paper, the MoG distribution has been considered to characterize the amplitude and 
the SNR statistics for wireless propagation. The parameters of the mixtures are evaluated by the 
means of the EM algorithm and the MSE and KL have been evaluated to challenge the proposed 
model’s accuracy. The proposed distribution enjoys both simplicity and accuracy and we have 
shown that the proposed pdf expression can accurately represent a wide range of both composite 
and non-composite channels. It should be highlighted that the adopted approach provides a 
generalized distribution for wireless communication systems where all channels can be modeled 
with the same analytical expression. Several analytical tools essential for the evaluation of 
performance analysis of digital communications were presented. This new model can be applied 
to various scenarios including, diversity systems, cooperative communications, and cognitive 
radio networks. 
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Appendix A 

MoG Parameters for Selected Scenarios 


The following tables provide the approximation parameters for all scenarios presented in the 
paper. 


Table I 

MoG Parameters for RL Fading Channel with £ = 3dB and C op t = 7 


i 

Wi 

fM 

Vi 

i 

0.24621 

0.76637 

0.1888 

2 

0.28164 

1.0954 

0.27482 

3 

0.15143 

1.5254 

0.39455 

4 

0.077823 

0.26911 

0.088815 

5 

0.025355 

2.0441 

0.59471 

6 

0.19662 

0.48573 

0.14066 

7 

0.020909 

0.11846 

0.050183 


Table II 

MoG Parameters for NL Fading Channel with m = 2, £ = 1 dB and C op t = 4 


i 

Wi 

Mi 

Vi 

i 

0.11009 

0.48576 

0.14483 

2 

0.14047 

1.3455 

0.34267 

3 

0.37935 

1.0845 

0.25099 

4 

0.37009 

0.77746 

0.19545 
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Table III 

MoG Parameters for NL Fading Channel with m = 4, £ = 1 dB and C opt = 4 


i 

Wi 

A H 

Vi 

i 

0.31126 

0.88351 

0.14927 

2 

0.13366 

1.2306 

0.25803 

3 

0.39008 

1.0756 

0.19479 

4 

0.165 

0.67541 

0.1405 


Table IV 

MoG Parameters for k — n Fading Channel with k = 1, n = 0.5 and C op t = 10 


z 

Wi 

(Jii 

Vi 

i 

Wi 

f-lii 

Vi 

1 

0.005263 

0.0052273 

0.0032608 

6 

0.1401 

0.39636 

0.11622 

2 

0.050502 

0.11458 

0.038823 

7 

0.013249 

0.020352 

0.0089878 

3 

0.18712 

1.4475 

0.3582 

8 

0.22111 

0.6682 

0.18341 

4 

0.21379 

1.0393 

0.24872 

9 

0.026951 

0.052855 

0.019739 

5 

0.054249 

1.8764 

0.48797 

10 

0.087675 

0.22199 

0.069679 


Table V 

MoG Parameters for k — p Fading Channel with k = 3, n = 1 and C op t = 5 


i 

Wi 


Vi 

i 

0.248 

1.2295 

0.29068 

2 

0.23197 

0.59067 

0.17727 

3 

0.017972 

0.24425 

0.099698 

4 

0.2658 

1.0976 

0.20736 

5 

0.23626 

0.86379 

0.16438 
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Table VI 

MoG Parameters for rj — ^ Fading Channel with 77 = 0.5, n = 0.2 and C opt = 14 


i 

Wi 

Mi 

Vi 

i 

Wi 

Mi 

Vi 

i 

0.092351 

0.21051 

0.061239 

8 

0.17858 

0.58555 

0.15984 

2 

0.030907 

0.034416 

0.012126 

9 

0.0052011 

0.00092585 

0.00068077 

3 

0.063814 

0.12078 

0.036005 

10 

0.0091546 

0.0047069 

0.0022762 

4 

0.092301 

1.8163 

0.41813 

11 

0.11756 

1.352 

0.28677 

5 

0.00030556 

4.5372 

0.47284 

12 

0.019634 

0.014034 

0.0055603 

6 

0.1304 

0.35591 

0.10026 

13 

0.029039 

2.5233 

0.61922 

7 

0.18817 

0.92188 

0.22362 

14 

0.042591 

0.067048 

0.021155 


Table VII 

MoG Parameters for 77 - n Fading Channel with 77 = 5, n = 10 and C op t = 3 


i 

Wi 

Mi 

Vi 

i 

0.32677 

1.063 

0.13602 

2 

0.32789 

0.89548 

0.099855 

3 

0.34533 

1.0116 

0.10407 


Table VIII 

MoG Parameters for k — 77 Shadowed Fading Channel with k = 1, p = 3, m = 3 and C op t = 4 


i 

Wi 

Mi 

Vi 

i 

0.40871 

0.86626 

0.18319 

2 

0.40942 

1.111 

0.22912 

3 

0.049081 

1.4037 

0.25915 

4 

0.13279 

0.61749 

0.15039 
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